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Abstract. We investigate the profile of random Polya trees of size n when only nodes of degree 
d are counted in each level. It is shown that, as in the case where all nodes contribute to the 
profile, the suitably normalized profile process converges weakly to a Brownian excursion local 
time. Moreover, we investigate the joint distribution of the number of nodes of degree d± and 
5-H ' d2 in the levels of the tree. 
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1. Introduction 



Consider the size of level k in a rooted tree, i.e., the number of nodes at distance k from the 
root. The sequence of level sizes of the tree is commonly called the profile of the tree. First 
■ investigations on the profile of random trees of given size (i.e., their number of nodes) seem to go 

back to Stepanov [SO]. A first distributional result as well as the relation to Brownian excursion 
local time was achieved by Kolchin [35] for the family trees of a Galton- Watson branching process 
conditioned on the total progeny. Different representations of the same result have been obtained 
later by several authors by approaching the problem from very different directions (random walks, 
queuing theory, general theory of stochastic processes, random trees), see [2"5I [501 15T1 I5"2"] . 

The relation between trees and diffusion processes has been studied by probabilistic methods 
as well. In this context, first and foremost the seminal papers of Aldous [TJ [5J [3] in which he 
developped the theory of continuum trees should be mentioned. Here, a metric space called 
C — ■ ' continuum tree is identified as the limit of several classes of random trees with respect to the 

Gromov-Hausdorff topology. This theory was further elaborated by Marckert and Miermont [37] 
and Haas and Miermont [29] . These approaches tell us a lot about the behaviour of large random 



trees. They imply limit theorems for global functionals. But they do neither provide limit theorems 
for local functionals like the profile of trees nor always convergence of moments. 

A different approach to the relation between trees and diffusion process was pursued in [42] . 
For a general overview on the relation between stochastic processes of combinatorial and therefore 
discrete origin and their continuous counterparts can be found in [43j . 

Not only the relation between trees and diffusion processes attracted much attention, but the 
processes themselves are of interest in their own right. A thorough overview on Brownian local 
times and related processes can be found in [46] . Explicit representations for the moments and the 
density of the one-dimensional projections of the local time of a Brownian excursion and related 
processes have been derived by Takacs [54] [53j [55] . Multi-dimensional analogues can be found in 
[2"71 12"8] . For results on density representations for related processes such as occupation times we 
refer to [T31 [33J E] • 

The profile of random trees has recently attracted the attention of numerous authors. A survey 
on the theory of random trees in general and the profile in particular can be found in [14] . Roughly 
speaking, the tree classes which have been studied can be divided into trees of height x logn (x 
meaning the order of magnitude) and trees of height x y/n where n is the size of the tree. Trees 
of logarithmic height are for instance binary search trees and recursive trees and variations or 
generalizations thereof. The profile of binary search trees and variations has been analyzed by 
Chauvin et al. [7] [5] and generalized by Drmota et al. [TO] and Schopp 08] . Hwang et al. studied 
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the profile of tree classes related to recursive trees, see [T5J [321 122] ■ The maximum of the profile, 
commonly called the width of a tree, in logarithmic trees was analyzed by Devroye and Hwang 
[12]. A related structure appearing in the theory of data structures are so-called tries. Their 
profile was examined in [TTJ [39l [41] . 

Trees of height x y/n are for instance trees obtained from conditioned Galton- Watson branching 
processes or Polya trees, i.e., rooted trees in the pure graph theoretical sense. The analysis of the 
profile of Galton- Watson trees was done in [IS]. The start of the profile, i.e., the behaviour close 
to the root, for trees as well as forests was investigated in [2H125]- Binary Polya trees have been 
studied by Broutin and Flajolet [5j [6J and general Polya trees in [17] . 

The joint distribution of two level sizes was addressed explicitely by van der Hofstad et al. [S5J 
and in [57] for Galton- Watson trees. Note that this question also arises implicitely when proving 
a functional limit theorem by showing convergence of the finite-dimensional projections as well as 
tightness of the profile, albeit the proofs of tightness utilizes the moments of the difference of the 
two level sizes. 

Recently, patterns in random trees were investigated as well. The questions considered here 
are the occurence of certain given trees as substructure of a large tree as well as the number of 
such occurences. The easiest pattern is the star graph. This amounts to counting the number of 
nodes of a given degree in random trees. First investigations in this directions were performed by 
Robinson and Schwenk [47] . In [16] it was shown for several tree classes (certain classes of Galton- 
Watson trees as well as Polya trees) that the number of nodes of given degree is asymptotically 
normally distributed, as long as the degree fixed. A phase transition occurs if the given degree 
grows with the tree size, see [35] for Galton- Watson trees and [33] for Polya trees. The analogues 
for general patterns instead of star graphs was carried out by Chyzak et al. [3]. 

Note that all results on patterns in trees tell us something about the number of occurences of 
a given pattern, but nothing about their location within the large tree. For Galton- Watson trees 
this question was settled by Drmota [13] for star graphs. The same question for Polya trees is 
addressed in this paper. 

Plan of the paper. In the next section we will recall some basic results on Polya trees and 
present the main results afterwards. The first result is Theorem Q] which states that the d-profilc 
(number of nodes of degree d at fixed distance from the root) view as a stochastic process weakly 
converges to Brownian excursion local time. To prove this theorem we will split it into two partial 
results, the convergence of the finite-dimensional projections of the profile (Theorem [2J as well as 
the tightness of the profile (Theorem [3j. To examine the joint behaviour of two different patterns 
we derive the covariance and the correlation (Proposition 12. ll and Theorem [4]). The results on 
covariance and correlation exhibit a surprisingly regular structure of the limiting object. We show 
that the correlation coefficient tends to 1 and derive the speed of convergence as well. 

The proof of these theorems will be carried out by means of generating functions. This will be 
described in Section [2] as well. This section ends with an introduction of the notation we will use 
throughout the paper. 

To proceed we need a singularity analysis of the generating functions together with a kind of 
transfer of the singular behaviour into the asymptotic behaviour of the coefficients in the sense 
of |20j . Section [3J provides some a priori estimates which are to be refined later. In Section [U 
we first present the proof of the one-dimensional analogue of Theorem [5] which is based on the 
refinement of the results in Section[3j The rest of the section is devoted to the refined analysis. The 
generalization to multiple dimensions is done in Section [5] In Section [6j we show tightness. This 
is usually a very technical matter (cf. the eight-page proof in |17j). Here we offer a considerably 
shorter proof by showing a more general result using Faa di Bruno's formula. The final section is 
devoted to the joint behaviour of the numbers of nodes of two different degrees within one level. 

2. Results and Notation 

2.1. Preliminaries. A Polya tree is an unlabelled rooted tree and thus it can be viewed as a root 
with a set of Polya trees attached to it. By the machinery of symbolic transfers described in [21] . 
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we easily obtain that the generating function y(x) of unlabelled rooted trees fulfils the functional 
equation 

V(x)=xexpfey^ , (1) 

a result going back to Polya [U] (cf also [55]) who also showed that y(x) has exactly one singularity 
p on the circle of convergence and that p « 0.3383219. Around its singularity, y(x) has the local 
expansion 

y{x) = 1 - b^/p - x + c(p - x) + d^/p - x 3 H , (2) 

with b « 2.6811266, as Otter [40] showed, and y(p) = 1. From the expansion, asymptotic estima- 
tions for y n can be derived by transfer lemmas (cf. |21|): 

Vn ~ ^-7=^— ( 3 ) 

2.2. Main results. We define by Ln (fc) the number of nodes of degree d at distance k from the 
root in a randomly chosen unlabelled rooted tree of size n. By linear interpolation, we create a 
continuous stochastic process L$(t) = ([t\+l- t)L [ n l) {[t\ ) + (t- [t\ )L^\[t\ +1), t > 

Theorem 1. Let 

and l{t) denote the local time of a standard Brownian excursion. Then l ( n d \t) converges weakly to 
the local time of a Brownian excursion, i.e., we have 

^ d) m>o^^-l( b At) , (4) 



2pb \2V2 



where C d = C + O (dp d ) with C = exp ^ 7 (^p- - l) J ra 7.7581604 . 



Remark. In [17] it is shown that the general profile of an unlabelled rooted random tree converges 
to Brownian excursion local time, i.e. 



t>o 



The normalising constant in Theorem [T] equals ^d^j^, where pan is asymptotically equal to the 
expected value of nodes of degree d in trees of size n, with pd — 1^7 /A see f° r example [3S] or 
[16] . To prove the above statement, weak convergence of the finite dimensional distributions and 
tightness have to be shown: 

Theorem 2. For any choice of fixed numbers t±, . . . , t m and for large d 



Remark. We will show this theorem by proving the convergence of the corresponding charac- 
teristic functions. It is well known (cf. []IJ|) that the characteristic function of C J^ h l f^f^) 

is 

r< n d r tJ^x exp ( — „ f b — x) 

m = 1 + ^E = 1 IV ^ — L dx (5) 

*JZ J -1 v^exp (^=) - iffi sinh (^=) 
where 7 is a contour going from +00 back to +00 while encircling the origin clockwise. 



4 



BERNHARD GITTENBERGER AND VERONIKA KRAUS 



A sequence of stochastic processes might not converge even if the sequence of their images 
with respect to every finite-dimensional projection does. Roughly speaking, in order to guarantee 
convergence in the sense of stochastic processes (i.e., when constructing a sequence by applying an 
arbitrary continuous bounded functional to the corresponding probability measures, this sequence 
must converge) the sample paths of the processes must not fluctuate too wildly. Tightness is 
a technical property of stochastic processes which guarantees this. The next theorem states a 
technical condition for the profile process which implies tightness (cf. [3] and [33] for the general 
theory). 

Theorem 3. There exists a constant c > such that all integers r, h, n the inequality 

E(L n (r)-L n (r + h)f <ch 2 n (6) 

holds. 

Remark. According to [4] Theorem 12.3] the inequality 

E \L n (r) - L n (r + h)\ a = O (/^(Vn) 0- ") 

implies tightness of the process l n (t) if a > and j3 > 1. In the theorem above we have a = 4 
and /3 = 2 and thus l n (t) is tight. We remark here that in [17j remark on p. 2050] the authors 
erroneously stated the bound O [(hy/n)^ . 

In order to examine the dependence of the numbers of nodes for two different degrees, say d\ and 
di (at the same level k), we will compute the covariance and the correlation. 

Proposition 2.1. The covariance Cov(X„ (fc)) of random variables Xn (k) and 

Xn (k) counting vertices of degrees d\ and g?2, with d\ ^ c?2 fixed, at level k — nn in a ran- 
dom Polya tree of size n is asymptotically given by 

Cov(A^)(fc), *«■)(*)) = C dl C d2 p d ^n (A [e~^ + - kV^) (l + O , 

(7) 

as n tends to infinity. 

Theorem 4. Let Xn dl \k) and Xn d2 \k) be the random variables counting the number of vertices 
of degree d\ and g?2> respectively, on a level k = Ky/n in a Polya tree of size n. Then the correlation 
coefficient is asymptotically equal to 1 as n tends to infinity. The speed of convergence is of order 
l/Vn. 

2.3. Description of the problem with generating functions. Proving Theorem [2 we will 
start with the one-dimensional case and then extend results to multiple dimensions. Therefore, 
we introduce generating functions y^{x, u), which represent trees where all nodes of degree d on 
level k are marked and counted by u. Note that we consider planted trees instead of 'ordinary' 
rooted trees, that is, we assume that the root node is adjacent to an additional node which is not 
counted. This assumption does not alter the tree structure, but allows us to treat the root vertex 
like a normal vertex, that is, a root of degree d has in-degree 1 and out-degree d— 1. Refining the 
decomposition of trees along their root, the yi (x } u) can be defined recursively: 

y[ f>(x, u) = y(x) + {u - l)xZ d ^ 1 {y{x),y(x 2 ), 2/(x d " 1 )) 
y$ 1 {x,u)=xe* P fc ^fi^ j, (8) 

where Z d {s\, s 2 , . . . , s d ) is the cycle index of the symmetric group 6 d on d elements, given by 

1 d 

where Xi is the number of cycles of length i in the permutation tt. 
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Examining two levels k and k + h simultaneously, we use the generating function y)jfl (x, Ui, u%) 
where all nodes of degree d on level k are marked by u\ and nodes of degree d on level k + h are 
marked by u-i- As before, x marks the total size of the tree. We get the recursive relation 

y i d) h (x,u 1 ,u 2 ) = y {d \x,u 2 ) + (tti - l)xZ d ^ 1 (y i h d \x,u 2 ), . . . , y i f ) {x d ~ 1 , uir 1 )) 

\i>i J 

In general, observing levels ki,k 2 = ki + hi,. . . , k rn — /c m _i + h m -i, we get: 

yO,h U -,hm-l • • • . U m) = Vhi,...h m -S X ' "2, ■ • ■ , U m ) 

+ (u- l)xZ d ^(yi d l hmi (x, u 2 ,..., u m ), y^l.u^ {x d -\ut\ • . • , u*- 1 )) 

Id) , \ / yk}ix, — ,h m -i i X% ' u l> ■ • ■ ' U m) 

y;>i 

These functions are related to the process (t) by 

F(Li d \k)=e u Ll?Hk + h 1 )=£ 2 ,...,Li d \k + Y,h l )=£ m ) 

[x n u^u^ ■ ■ ■ u l ^\y^ hu ^ hm _ x (x, Ml, ... , u m ) 
[x n ]y(x) 

For the computation of the covariance of the numbers of nodes of degrees d\ and d 2 we will utilize 
the functions 

y^' d2 \x,u,v) = y(x) + (« - l)xZ dl ^(y(x),y(x 2 ), . . .^(x^ 1 )) 



+ (v- l)xZ d2 ^(y(x), y{x 2 ), y^ 2 " 1 )) 



2.4. Notations. For proving the theorems we will carry out a singularity analysis of the generating 
functions. All generating functions are in some sense close to the tree function y(x) from (Q}. 
Therefore we will use the differences to y(x) and related functions: Set 

w[ d \x,u) = y[ d \x,u) - y[x) wf 1 '^ (x, u, v) = y^ uda \x,u,v) - y(x) 

E W> (x,u) = J2 ^ ( f X) 4* M (w) = E ^ 1,rf2)( f' M>I) 

i>2 i>2 



7fc J (^w) = — y y k >{x,u) 

We further introduce some domains, depicted in Figure [T] 

A= A(r),8) = {zeC||z| < p + 77, | arg(z - p)| > 9}, (12) 

A e = A c (6) = {z e C\\z - p\ <e,|arg(z-p)| > 9}, (13) 

9 = 6(77) = {zeC||z| <p + ?7,| arg(z - p)\ ± 0}, (14) 

S fc =Z k (rj) = {ze C\\v\ < l,fc|«-l| < fj}, (15) 

with e, 77, fj > and < < §. 

In all the proofs in the subsequent sections we will assume (even without explicitely mentioning) 
that 77, 9, e are sufficiently small for all arguments to be valid. 



(i 
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Figure 1. The regions used for the proofs 

3. The local behaviour of y[ d \x,u) and y < f dl ' d2 \x,u,v) - A priori bounds 

In order to analyze the local behaviour of the generating functions we will first derive a priori 
estimates for w^, u>ji, dl ' d2 ' ) and the related functions which will be used frequently in the sequel 
to derive the needed refinements. 

Lemma 3.1. Let |x| < p 2 + s for sufficiently small e and \u\ < 1. Then there exists a constant L 
with < L < 1 and a positive constant D such that 

\w[ d) (x,u)\ <D\u-l\-\x\ d -L k 

Proof. We will only provide a short sketch, since the proof is similar to that of [T71 Lemma 2] . 
For k = we have 

\w (d) \(x, u)\ = \u-l\- \x\ ■ |Z d _i(y(x), y(x 2 ), . . .)| < \u - 1| • D ■ \x\ d 

0{\y{x)\^)=0{\x\^) 

The result for general vJi (x, u) follows by induction. Starting with the recurrence relation 



Wfc+ifou) = y{x) exp w (d \x,u) + k ■ 



we use the trivial estimate \wk{x, u)\ < 2y(\x\) which is valid for |x| < p and \u\ < 1, the convexity 
of v( x )/x on the positive reals, and some elementary estimates for e x . For the precise details see 

021- ' □ 

Corollary 3.2. For \u\ < 1 and \x\ < p + e there is a positive constant C such that (for all 
k>0,d>l) 

\El d \x,u)\<C\u-l\L k . 
Proof. Same as proof of Corollary 1 in [17] . □ 
Corollary 3.3. Let u <E and x € 0. Then 

Y,li d \x l ,u>) = 0(L k ). 

i>2 

Proof. As i > 2 the functions 7^ (x l , it 1 ) are analytic in the whole region and T^ d \x,u) := 
T, i > 2 ~fi d) (x\u t ) = T, n ,mVnik xn y m with positive coefficients 2/^ m , we have \T {d \x,u)\ < T {d) (\x\, \u\) 
where the right-hand side is monotone in |x| and 

Now let x > and < u < 1. Using Taylor's theorem we get 

4 d \x, «) = (« - i)ri d) (x, 1 + 0( U - 1)) > ri d \x, u). 

In view of Corollary 13.21 this implies for all x g and it € 3* the estimate |rj^(x, u)| < 
r^dx], < CL fc for some positive constant L < 1. □ 
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In a similar fashion we obtain the analogous results for w^ 1,d ^: 

Lemma 3.4. If \x\ < p 2 + e for sufficiently small e, \u\ < 1, and \v\ < 1, then there exists a 
constant L with < L < 1 and a positive constant D such that 

\w[ dl ' d2 \x lU ,v)\ < DL k {\u~l\ ■ \x\ dl \x\ d2 ) 

Corollary 3.5. For \u\ < 1, \v\ < 1 and \x\ < p + e there is a positive constant C such that (for 
all k > 0,d > I) 

\x{ dl ' d2 \x,u 7 v)\<C(\u-l\ + \v-l\)L k . 
Corollary 3.6. Let u G E/., v £ E/., and x E 0. Then 

Y,l[ dl42 \x\u\v*)=0{L k ). (16) 

i>2 

4. The one dimensional case 

Our main goal is to prove the following theorem from where the main result follows by integra- 
tion. 

it 

Theorem 5. Let x = p(l + — ), u — e^™ , k — \Kyfn\ and d be a fixed integer. Moreover, assume 
that |args| > $ > and, as n — > oo, we have s = 0(log 2 n), whereas k and t are fixed. Then, 
wif* (x, u) admits the local representation 

{d) . A C d p d it^Se-^ b ^~ s 



wi '(x, u) — — =- ■ — = — d — — (17) 

The one-dimensional limiting distribution. Let us first assume that Theorem [5] holds. Then, 
to prove Theorem [2] in one dimension, we need to determine the characteristic function 



— U*i*,e*)^ (18) 



where the contour r = 7 U T' consists of the line 

7 = {x = p(l - 1 + tT )\ - D log 2 n < t < L>log 2 n} 

n 

with an arbitrarily chosen constant D > and T' is a circular arc centered at the origin and 
closing the curve. The contribution of V is exponentially small since for x € V we have 
-^-\x~( n+1 ^ \0{nie~ log2 ™) on the one hand whereas on the other hand \y^ d \x, e^™ )| is bounded. 
If x € 7 the local expansion (|17p is valid and thus, inserting into (fT5|) leads to: 



fc(«0 , 1 



lim (pi ' (t) = lim 

n— foo ' n— foo 



27riy„ 



(d ) , cfe /" , , dx 



-2-Kiy n 



_ log n , /— — 

C d p d+n nV2 f t^/^e^"-^^) 1 



n— »-oo 



1 + lim — Lir , , — / ■ — v ■ e s ds 

bj2pH J e (^S) _ s - mH ^E) p n n 



1— ilog 2 n " VP b 

where ip(t) is the function given by (0. 

Now let us turn back to the proof of Theorem [5j 
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4.1. The local behaviour of yr (x,u) - refined analysis. Now we will refine the a priori 

estimates of the previous section. First we show that the first derivate jjf\x, 1) is almost a power 
of y(x). Afterwars we will derive estimates for the second derivative and then obtain a power-like 
representation for lol (x, u). Finally, utilizing the recurrence relation for «;[ (x, u) we will arrive 
at the desired result (fT7|) . 

Lemma 4.1. For x 6 <3 (where ?y > is sufficiently small) the functions Jjf\x) can be represented 
as 

li d \*) :=7i d) (x,l) = Ci d Hx) y (x) k + d , 

where the functions C% (x) are analytic and converye uniformly to an analytic limit function 
(x) (for x £ Q) with convergence rate 

Cf{x) = C {d \x) + 0{L k ) 

for some < L < 1, and further C^(p) = Cdp d , where Cd is the constant given in ((4]). 

Proof. We define the functions cf {x) := 

We prove the analyticity of the functions 7^ {x) by induction: 

^(x) = xZ d ^{y{x\ y(x 2 ), . . . , y^ 1 )) = xO^xf- 1 ) 
is analytic in by previous arguments, and so is 

_ xOiyjx^- 1 

y(x) 

The step of induction works like in [T7] , as the 7^ fulfill the same recursion as the 7fc 



C o W= 1^ J =0(l) (19) 



7*Ji(a,u) = it***- 1 

E£^( 



<9?r 



= ^ 1 (x, M )^ 7 f(^,uV- 1 , (20) 
i>i 

and for u = 1 

7SiW=y(^)7i rf) W+j,(x)rf (x), 

with r^^x) = llf^ix 1 ), which is analytic for \x\ < ^fp and hence in 0. Applying the induction 

i>2 

hypothesis, this proves the analyticity of 7^ (x). Solving the recurrence, we obtain 

k-l 

% 



^ = y (xM\x) + J2y^) k -^ ( f ) ^) 
1=0 



and hence, the analyticity of 7^ implies the analyticity of the functions C k (x) in 



We now have to show that the functions (C^ (2;))fc>o have a uniform limit C^ d \x) but this works 
analogously to |17[ Lemma 3]. 
Finally, note that 

fe>0 fc>0 

where dn is the total number of vertices of degree d in all trees of size n, and D^(x) is the 
corresponding generating function, introduced in e.g. [47] . On the other hand, 

E 7fK 1) = E( G(d) (-) + 0(L k ))y(x) k = C ^[ X lf x f + 0(1) 

fc>0 fc>0 ^ ' 
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and therefore 



We know that 



y(x) 



x—>-p 



1 - y(x) 



(cf. [47l Eq. (36)] or [36]). Schwenk [49l Lemma 4.1] computed the limit of the cycle index in the 
numerator. In his proof he provides the speed of convergence as well. In fact, [49J Eq. (32)] says 
that 



Z d 



y(x) y(x d ) 



ex p E - 



1 / y(x) 



- 1 



1 \ x 



with A = supo^ i - l) = Thus xZ d ^(y(x), . . . ,y{x d ^)) = x d F(x) + O (dr 2d+1 ). 

Note further that D^ d \x) = O since there are no nodes of degree d in trees of size less than 

d+ 1. This implies C^(p) = C d p d with C d = C + O (dp d ) and C* as in Theorem [Q □ 

Lemma 4.2. There exist constants e (d) , 6^, rj( d) > and 6^ < § smc/i i/ioi 

| 7 i d) (x, U )| = 0(|y(x)r d ) 

uniformly for x £ A £ and afHfc. 
Proof. For I < k we set 



- sup 



7£ 



First we derive the following inequality, using the recurrence for y k 



(<*). 



cxp ( E -yi^fcW) - y 
\i>i * 

Or) exp (j^jwfix^u^ 
<|y(z)|exp \\w^(x,u)\+J2^e d H 



< I y (a;) I exp 



| 7 f \z, l+d(u- 1))||« - 1| ^7 ^"V, 1 + " l))l 



V 



<C< d) |u-l| 



i>2 



J 



with < d < 1 and thus 1 + {)(u l — 1) <G 3&. To get an estimate for the second term, we use that 

|u* - 1| = |1 + u H h - 1| < i\u - 1| as |u| < 1 and hence < |u - 1| < 2. Further 

we use \yf \x\ 1 + - 1))| < \lf\x\ 1)|, |j/(x)| < 1 and Corollary E3] to obtain 



lyffxMI < |y(^)|exp (cfV-i| + o(^)) 
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Using recurrence (j2"0)) leads to: 



cffi = SIl P 

u£3 k 



Vl^ix, u) 



y{x) 



2/(a0 



< e c?r^+o(^) ( cw+o(^)) 

= C i (d) e^ > ^(l + 0(^)) I 



(21) 



where we used Lemma liTTl to get | J2i>2 l^f\x l ,u l )u k 1 \ = 0(^2 i>2 \y{x l )\ e,+d ) and hence 



sup 

xt£A c 

We now set 



Note that 



( d )/^i „,i\„,i—l 



y(xY 



( y{* 



i>2 



V y( x ) 



= 0(L e+d ) = 0(L l ). 



j>0 



7o d) (^ u )i _ ,x z d-i(y(x),y(x 2 ), . . .y{x) d *) 



y{x) 



y{x) 



1 = 10(^)1 = 0(1), 



hence = sup | 7 °^)d" ) I = 0(1), too. Thus we can choose rj > such that e 20 ^'? < 2. For 
fixed fc we get: 

Q{d) < ^d) -q (1 + o^^gae^ooci < 2 (7 oCo = 0(i). 

i<l 

The second estimate is clear by the choice of r\ and by I < k. The first inequality can be obtained 
from (f2Tj) by induction: 



C[ d) < C { d) (l + 0(L ))e°o d> i < C { Q d) H(l + 0(L>))e 2C o 



Cg> 1 <CV>ei c ?\l + 0{L 1 )) 



3<1 



= JJ(1 + 0(L j ))(l + 0{L £ ))C { d) e 2C o d)c oni ey]l 

3<l + l 



<2fc<2 



For the second derivatives with respect to it of yi (x, u), 7! (x, u), we find that 
Lemma 4.3. Suppose that \x\ < p ~ rj for some r\ > and \u\ < 1. XTien 

7 P ] M = 0(y(\x\) k+d ) 
uniformly. There also exist constants e, 8, f) such that uniformly for ugHi tmd i£ A e 

7 f [2l (x, W ) = 0(fcy(kl) fc+d ) 



□ 



(22) 



(23) 
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Proof. Derivation of (|2"U)) leads to the recurrence 



with initial condition 7g d ^ 2 '(x,u) = 0. 

For |jc| < p — r), for some r\ > and for |u| < 1 we have (7^ (£!m)| < 7^ 1)- Thus, in 

this case we can restrict ourselves to non-negative real x < p — r/. 

By using the bounds 7^ (x, 1) < C^y(x) fc+d from LemmaOJ Ei^fc^W) = °( Lk ) from 

Corollary [33] and the induction hypothesis J^^ 2 \x, 1) < D k d ^ y{x) k+d , we can derive the following 
upper bound from the above: 



r ,i\k+d\ 



£e<r^) fc+d ) 



y;>i / z>i i>2 

< y(x) (cfy{x) k+d + 0(L fc )) 2 + (£ fc( 

< y(x) fc + d+1 [cf )2 y(a;) fc+d + cf >0(L*) 



//(^) fc+d i , c y2(i-i) y ( xi)k+d 



i>2 



y{x) 



k+d 



i>2 



y{x) 



k+d 



< y(x) k+d+1 {C[ d)2 y{p - rj) k+d + D k (l + 0(L k )) + 0{L k )). 



Consequently 



Dfl = (C { k d)2 y{p + D k (l + 0(L k )) + 0(L k )) 



\k+d 



which leads to — 0(1) as k — >• 00. 

To prove the second property we use the same constants e, 9, 77 as in Lemma 14.21 and set: 



{dm, \ 
D e = SU P 7 Mt+d ' 



We use the already known bound \"/p d \x,u)\ < C^\y(x) k+d \ and by similar considerations as 
in the proof of Lemma |4~21 we get: 



12 



BERNHARD GITTENBERGER AND VERONIKA KRAUS 



£7+1 = SU P 

u£E k 



J2 7£ d) (x I ,W )// 

i>l 



i-1 



i>l i>2 



y{x) 



< D^e^Hl + 0(L e )) + C™e eW i + 0(L 1 )) 



with a 



(rf) 



S c(d) *(l + 0(Z/)) and = C'^e^'t + 0(L e )). Thus 



^fe < afe-i("fc- 2 (- • • \ a o Uo + Po ) ■ ■ ■)P k -2) +Pk-1 
k—1 k—1 

TT ^ d ) _l „,( d ) r>( d ) 



=E/f n 

3=0 i=j+l 

< fcmax/?f ) e C ' (d)c ]^[(l + 0( J L 2 )) 



i>0 



= o(*), 

which completes the proof of the Lemma. 
Lemma 4.4. Let e,9,fj and cjf\x) be as in Lemma \J^.l\ and Lemma \4-2\ Then 



w {d \x, u) = C {d \x){u - l)y{x) k+d {l + 0(k\u - 1|)) 



uniformly for x G A £ and u € Sfc. Furthermore we have for \x\ < p + t] and \u\ < 1 
4 d) (z, u) = C {d \x){u - l)y(x 2 ) k+d + <D(\u l\ 2 y(\x\ 2 ) k+d ), 



where the analytic functions c£ (a;) are given by 



□ 

(24) 
(25) 



i>2 



y(x 2 ) 



and /iave a uniform limit C^ d \x) with convergence rate 

c\ d \x) = C {d \x)+0{L k ) 

for some constant L with < L < 1. 

Proof. To prove the first statement, we expand lol (2;, u) into a Taylor polynomial of degree 2 
around u = 1 and apply Lemmas 14.11 and 14.31 

To prove the second statement, we again use Taylor series. Note that for i > 2 we have 
\x % \ < p — r\ if \x\ < p + 77 and ?y is sufficiently small. We get 

<4 d VX) = Cf V)(^ - lM^) fe+d + 0(\u l - \\ 2 y(\x l \) k+d ) 
and consequently 

T<g>{x,u) = £ - l)y(x l ) k+d +0(\u~\\ 2 y{\x 2 \) k+d ) 

i>2 

= (« - l)cf > (zMz 2 )^ + - l| 2 y(M) fc+d ), 
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where we used the property that 

it* - 1 y(x l ) k+d 



i>2 



i(u — 1) y(x 2 ) k+d 



(1+tH hit 1 " 1 ) y{x l ) k+d 



i>2 

Cf\x)+0{cf\x){u-l)) 



V 



represents an analytic function in x and u, and thus its leading term, as u — > oo, is our function 
C^fx). Finally since Cf\x) = (x) + 0(L k ) it follows that djf \x) has a limit ^(i) with 
the same order of convergence. □ 

Lemma 4.5. For x £ A e and meSj (iirci/i t/ie constants e, 0, 77 as in Lemma \4-.2ty we have 

<<t),. A _ (u-l)y(x) k + d Cf(x) 



w k '(x,u) 



Proof. wjp(x,u) satisfy the recursive relation 



(d) 



(x,u) = xexp [ £ -y k d) (x\u l ) I - y(ac) 



y(x) 



y(x) (exp (wj^ (x, u) + (a;, tt)) - l) 



and further, since by LemmaHUit follows that zZ k d) (x,u) = 0(w^ (x, w)i fe ) = 0(^ d) (a;,u)) (for 
brevity, we omit the variables now), 



(d) 

<+i = 1/ 



+ 4 d) ) + K + 2 Sfc } + o ((«,«> + 4 d) ) 3 ) 



y(«|W> + E«>) ( 1 + 



2 

.,(<*) 



+ 0((«;W+sW) a 



o(4 rf) ) + o(( w i d) ) 2 ) 



From there, we obtain 

2/ 

(d) 
<+i 



.(d)- 



(<') 



^ (i + ^ + 0(4 d) ) + 0((4 d) ) 2 ) 
^fi-^+^(4 d) ) + ^ d)2 )) 



,.(<*) 



,(dr 
„( rf ) 



This leads us to a recursion 



,fc+i 



(d) 

4+1 



„(<*) 



(d) (d) 
w fc w fc+i 



0(w k d) y k ) 
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which we can solve to 



(d) (d) 
w k w o 



fc-1 v (d) 

E 

^=0 



fc-1 



(d) (d) 

W £ W £+l £=0 
fc-1 „(d) e + x 



E 



fc-1 v (d) * 



■0(X>JV) 

f=0 



(d) 



i-< ] E 



(d) (d) 

w £ «>/+! 



(26) 



where we used that 



0(1/) and that by Lemma WMw [d) = 0(yw^\) = O(y k w { a) ). Again 



(d) 



,fc„„( d h 



we apply Lemma l4~4l and (jSJ to obtain 



.,(<*) 



fc-1 v (d) 



E 

1=0 



^ - %(x 2 )^ + C% - l| 2 2/(|x| 2 )^) /+1 

= (u— ±)xZd-i > — — tt; V 

^ cfc\%yW+^{u mi + 0(£\u 1|)) 



sZ,,-! ^ cf Vx 2 /+ rf + C% - l| 2 y(M 2 )' +d ) 
V(x) d S cf ^y^^l + 0(£\u - 1|)) 



y{x) 



■fc-i 



E 



y(x 2 ) 



2\i+d 



fc-1 

E 



0(| M -l| 2 j/(|x| 2 ) £+d ) 

J fi{d) r (d) v ( x )l+d 1 r (d) r (d) , Y+d 

W W+l yV ; £=0 4> W+l2/W 



j/ (a; 



^gd-i 
y(^) d 



ci d) +0(| M -l|), 



— C] rf) j/fc 2 )** 1 ^ 0(\u - l| 2 y(N 2 )^ +d ) 



E 



^ ^(^(O y(x) e+d ^ ni d )ri{ d ) ,,(^\l+d 

°£ W+l yV 7 fcO W L^+li/W 



=0(|«-1| 2 L«) 



1 + 



(l + 0(*|u-l|)) 



(d) 



y(a; 2 ) 



fc-i 

E r ,!,/, r ,i./! ,,( ,.)'+•' 

f=0 °£ °£+l yv y 



2\l+d k-} £<( rf ) 



E 



2/0 



„2\l+d 



r (d) r (d) v ( x y+d 

=0 °£ ^£+1 

=0(L«) 



0(%-l|) + 0(| u -l| 2 ) 



where ci denotes the first sum. Note that X y* x y = 0(1). 



Now turn back to ([26]) and observe that w ( d) = £>(% - l|y(a;) d ) = 0(y(a:) d ) = 0(1) if 

fc|it — 1| < fj. Thus, we obtain the following representation for w' yd \x, u): 

(d) fc 

„..(<*) _ 



l-4 d) (x)-^i^ + 0(| U -l|) 



(d) fc 

(d) 

2 l-j/ 



We use the expressions 



°fc+i - 



(d)2/0 



„i\fc+d 



i>l 



y(a; 



,fc+d 



and(7f=^^ 



i^ifc+d 



j>2 



(d) 2/(^) 

y( X 2)k+d' 



which are consequences of Lemmas 14. 1[ Equation (|2T)|) and Lemma 14.41 to obtain 



ci d \x) = (Cl%(x)-Ci d H X ))l y{x) 



y(x 2 ) 



fc+d 



(27) 
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This provides the telescope sum: 



fe-i n (d) r (d) 



v ( x )d r ,(d) r ,(d) 

yy > e=o u e W+i 
xZ d -i I 1 1 

y(x) d \c (d) c (d \ 



(29) 
(30) 



and hence, since C^ d) = -$^r = ^f^f, we get 



y{x) d y(x) d ' 

1 _ (<i) r \ _ xZd-Ay^ ■ ■ ■,y{x d - 1 )) 
~ Ck {X) y(x) d ci d \ x ) ' 

which yields the result. □ 



It is now easy to proof Theorem [5j With x = p(l + fr), u = e^, d and £ 7^ fixed, k = K^/n 
and ([2]) we obtain the expansions: 



it 

u — 1 ~ — = 



l-l/(x)~6 V ' ' PS 



y(a;) fc ~l-fc&W-^ + e 



Since the functions cl (a;) are continuous and uniformly convergent to C^ d \x), they are also 



uniformly continuous and thus c[ d \x) ~ C^(p) — Cdp d - This leads to 



( d )i \ 
w\ '{x,u) 



it s i n d„ — Kb 
7^ C dP e 



1 l-e-^v 

2 >,. /3 



-sitCdp d e 



c dP d 



itC d p d I 
2b^p ' 

it 



1 



-Kb\/—ps \ 



\K,b^=p~i _ itC dP d 



^sinh(i«V^) 



(31) 



5. Finite dimensional limiting distributions 

First we consider the case m — 2. The computation of the 2-dimensional limiting distribution 
shows the general method of the proof. Iterative applications of the arguments will eventually 
prove Theorem [2J 

it 1 it 2 

Theorem 6. Let x = p(l + ui = e^,^ = e^™, k= n^/n and h = r/y/n. Moreover, assume 
that I args| > O > and, as n — » 00, we have s = C(log 2 n), whereas k, t\ and t2 are fixed. Then, 
for large d, w^^XjU) admits the local representation 
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wl' h {x,u,v) 



c dP d 



ito 



iti^/ — se 



= >-±i^sinh(i«V^)) 



/^ e (-h£ b V-P s ) 



- 1 / [— -i-rebV— ps) 

■lt\y/ — Se y 2 v F J 



(32) 



e( i „ t ^, _ u 1 c^_ s . nh ( 7 KbV=m) 



sinh(±£6 v / =ps)) 



Proof. Note that y^ h {x, iti, 1) = y^f\x, u{) and V/f^ix, 1, u 2 ) = yk+h(x, it 2 ). Considering the first 
derivative, we denote by 

d (d) , \ (d)M/ \ 

Q^yk,h( x ' u ii u v = : %,h {x,ui,u 2 ), 



and by simple induction, we observe that: 



7 i5 [tiii (x, 1, i) = 7 nx, i) = 7rw - crww^) 

7 i5 [ - ] (x, 1, 1) = 7 W fc(jB> i) = 7 W h(a;) = C^C*)^ 



(33) 



As l7[. d l [ "* 1 ( a; : u i> w 2)| < j^ Ui \x,l,l) for i = 1,2; m G S fc ,u 2 G Efc+ft) and \x\ < p it follows that 

|7^ [tll] (x,tii,u 2 )| = C(y(x) fe + d ) and 7 ^ [ " 2] (x, u x , tt 2 ) < 0{y{x) k+h+d ) in the same regions. To 
be more precise, we can prove the following analogue to Lemma [ 



Lemma 5.1. There exist constants e, i?, 771 , 772, such that or x € A e ,m € Sfc and u 2 € Sfc+/, 
7 g M (*, u x , « a ) + 7 <5 M (z, u x , « 2 ) = 0(\y(x)\ k+d ) 

Proof. Set 



e.h 



C, 



(d)[u 2 ] 

e,h 



sup 

u 1 £E k ,u 2 £5 k+h 



sup 




As in the proof of Lemma T4.3[ we apply Taylor's theorem (in two variables) to get 



m+i,h( x > v 'i> u 2)\ = exp \w\ d) {x,u 1 ,u 2 )\ 



\w ( e d) (x\u\,u 2 )\ 



i>2 



< \y(x)\ exp (ifl [ui] {x, 1 + ft (ui - 1), 1 + d 2 {u 2 - l))(ui - 1) 



+ 7 <f Ua] (or, 1 + ft (ui - 1), 1 + M^2 - l))(ua - 1) 



E 1 + - !). 1 + M< 1)) 



K - 1) 



i>2 



+ ^f U2 ' {x\ 1 + 0i K - 1), 1 + Mu l 2 - 1)) 



u\ - 1 



(d) 

Vt+i 



U(^«i.«a)| < l^)l ex P (cS )M k - %W +<i + Cg [u2l |« 2 ~ %(^ +d + 0(^)) 
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where we use that, for i > 2, 

| 7 g M (^l + ^iK-l),l + ^(4-l))l < bg^V, 1,1)1 and 

| 7 g [ " 2 V. i + M< i), i + M4 - i))l < l7i?" 2 V, 1, 1)1- 

By using recursion @ and Lemma |4.2[ we obtain 



sup 



< exp (C, 



, r (d)[u 2 ]V2 

''" k e ' h k 



(d)K], 



(<*)Mm 



(<i)[ui] 



+ 0(L e 



k+d 



= exp (cg^f + (1 + 0(L<)) , 

and analogously 

. (d)[ ui l .-(JJIuj] 

We choose 771 and r\ 2 such that e ° y °- h ni ^o,h T n < 2. Then, by induction we get 

C& )IUl] < <2 [ui1 + 0(L J '))e 2c -(°or i] '»x+^ [ " ! '^)i < 2C$ M cd = 0(1), 



< c gK] + 0( ^ ))e 2 Co( c-'-S 1+ c-'-S 2 )4 < 2d d) h ^ ] c = 0(1). 

Note therefore that 

(J)W . „„„ a;^d-i(2/i d) (a;,W2),...,2/l rf) (a; d ~ 1 ,u2 _:L )) 



0,/i 



Ci 



MM 



0,h 



sup 



sup 



y{x) d 



0(1) 



j^ d) h [u2 \x,u 1 ,u 2 ) + (ui - l)x-^Z d _ 1 (y h (x,u 2 ), 



= 



sup 

i6A e 
1 uiSH fc ,ii2GHfc + fe 



y(x 



,h+d 




a 



Let 



(d)[2u 1 ] l 
<k,h 



d 2 



%,h ix,u 1 ,u 2 ) := -^y kyh {x, Ul ,U2) 

1 " 1 - ii-i 1 , @ „ (d) / \ 

Tk,h {x,ui,u 2 ) := Q^y kih {x,ui,u 2 ) 
d 2 



(d)[u!U 2 ], n {d) , , 

j kh {x,u 1 ,u 2 ) := duidu yk,hKX,Ui,u 2 ) 

7 £f ] (*,ui^ 2 ) := 7 ff Ul k^"i,"2)+7i d r 2 k^«i I " 2 )+7S [tllH (x, U i W2 ) 
Lemma 5.2. For \u\\ < 1, | U2 1 < 1 a^d for \x\ < p — r\ for some r\ > 

1 if 2 ^(x,u 1 ,u 2 ) = 0(y(\x\) k + d ) 
7 g^(a ; ,u 1 , U2 ) = 0( y (|x|) fc + ft + M - 1 ) 
7 S [2tl2] (x, Ul , U2 ) = 0(j/(|x|) fc+ ' l+d ) 
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uniformly. Furthermore, for i 6 A e ,«i G sj. and u 2 G S 



k+d\ 



Proof. The proof of the first statement is identical to the one of Lemma 14.31 as we can derive 
identical recursive relations for ^/jfj^ 2ui \x,Ui,U2) and T^' 2 " 2 ' (x, ux, u 2 ) and a similar one for 

{d)[u 1 u 2 \ / \ 



i>l / \i>l 

(d)[«i«2]^i „,j ^^(t-l^t-X) 



i>l 



We then prove the statement inductively with the following initial conditions (note therefore that 
£-Z n (si,...,s n ) = i Z n -i(sx, . . . , s n -i) (cf. [361 Chapter 2, p. 25]): 



T^^^ux,^) < = 0(y(x) h+d ), 

d_ 
du 2 



^f^\x,u 1 ,u 2 )=x^Z d _ 1 (yi d \x,u 2 ),...^\x d -\ut 1 )) 



d-1 



8 



^-^d_i(si, . . . , sd-i) 

r—l r 
d-1 ^ 

^ -^d-r-l( s li • ■ • i s d-l-r) 



id)/ r r \ r—l 



s i =y h (x i ,u' 2 ) 



(d) / r r\ r—l 

l h {x ,u 2 )ru 2 



r—l Si=y h {x\u\) 

= 0{y{x) h+2d - 2 ). 



For the proof of the second statement we define for I < k 



D 



(d) 
t,h 



sup 



(d)[2]/ \ 



y{xf 



as in the proof of the second part of Lemma 14.31 We use the estimate 



y^ h (x,u u u 2 )\ < | y(a; )|exp(c2^| + Cg M | + ^)) , (34) 



which we obtained in the proof of Lemma 15.11 

From the recursive description, we can derive the following by applying known bounds from 
Lemma 15.11 and from the previous statement, and from (|3"3"f and (|3~i]) . similar to the proof of 
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Lemma 



D 



e+i,h 



sup 

xGA E 
uiGH fc ,u 2 GH fc+fe 



sup 



li+ih (^,Wi,W2) +7^ 1 [ ^ t2j (x,ui,u 2 ) +7)+il 1 " 2l ( a; '' Ul '" 2 ) 



U)[2u 2 ] 



(d) [u 1 u 2 ] , 



ife+d+l 



E^i(E i >i7g Kj (^%ui,4)4- 1 ) 2 + re=i(E l >i7^ lUr H^%^,4> 



2/(3) 



(d)K 



,t \„,i— 2 



exp 



fC 1 T + C S W T + °( if ) 
x ((Cg [lll] ) 2 y(a;) fc + d + C\f u ^C$}^y(xf +h + d + (C$ [u2] ) 2 y(x) k+2h+d + £>g> + 0{L l ) 

+ exp (cg M f + Cg [U2] f ) ((< )M ) 2 + Cg^^i?" 21 + (Cg [u21 ) 2 + 0{L*))) 



As in the proof of Lemma T4.3I we get 



fc-i 



fe-i 



(d)> 



i=j+i 



0{k) + 0(D^>) 



(35) 
(36) 



It remains to prove that -Dg ^ = C(M 



7$ laui] 0r,ui,« 2 )=O 



7$ [2ua] (*,«i,«a) = 7 i d)[2] (x,u 2 ) + ^^(y^fou.,), . ...yW^-i.u-- 1 )) 



1=1 j=l 

+E^-i(^ ) )7r t2] (^,«')^- 1 

+Ez d _ ; _ 1 (,f)a-i)7i d) (x',^) U '- 2 
;=i 

0(%(x)' l + d ) + 0(y(a;)' l + 2d - 3 + OQiy{x) h+2d - 2 + 0{y{x) h+2d - 2 ) 
0{hy{x) h+d ) 
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and 



^ f u ^ ] (x, Ul) u 2 ) = A-Z d ^(y (d \x, u 2 ), . . . , yi d \x d -\ ut 1 )) = 0(y(x) h + 2d - 2 ) 



D (d) _ 
u o,h — 



sup 

x£A t 

u 1 es k ,u 2 eE k+h 



(d)[2ui], s , (d)[2u 2 ]z \ , (d)lu 1 u 2 ] r \ 

Toh (x,u 1 ,u 2 )+%,h ( x , u i, u 2)+r h > {x,u 1 ,u 2 ) 



y{x) a 



< sup 0(hy(x) h + y{x) h+2d - 2 ) = 0{h) 
u 1 es k ,u 2 es k+h 



□ 



Remark. Obviously, for x € A e and «i £ ^,^2 £ Sfc+fc the statement also holds for the partial 
derivatives: 



l { k T Ul] {x,u u u 2 ) = 0{ky{x) k+d ) 

iiT U2] ^ u ^ = °« k + h)y(x) k+h+d ) 

iiT lU2] ^ Ui < M2 ) = + h )y( x ) k+d ) 



Lemma 5.3. For x e A e ,Ui € and u 2 € ^k+h, with the same constants as in the previous 
lemmata, we can approximate 

w {d l{x,u u u 2 ) =C {d \x){ Ul l)y{x) k+d + C {d l h {x){u 2 l)y{x) k+h+d 
+ 0((k + h)y(x) k+d (\ Ul - 1| 2 + \u 2 - 1| 2 )). 

Furthermore 

XW(X,U U U2) =ci d \x 2 )( Ul - l)y(x 2 ) k + d + C^ h (x)(u 2 - l)y{x 2 ) k+h+d 
+ 0(y(\x\ 2 ) k \ Ul -l\ 2 + y(\x\ 2 ) k+h \u 2 -l\ 2 ) 

Proof. For the first statement, we expand w^ h (x,ui,u 2 ) into a Taylor polynomial of degree 2 
around u\ = u 2 = 1 and obtain 

Wfc^wi,^) = 7i d) (^)(«i - 1) +l ( k+ h (x)(u 2 -1)+R 

with \R\ < 1 (tS' 2 " 11 ^, 1 + - 1), 1 + tf 2 («2 - l))(«i - I) 2 

+ 27S [U1U21 (^, 1 + Mui - 1), 1 + ^ 2 («2 - l))(«i - 1)(U2 - 1) 

+7S [2 " 2l (^, 1 + - 1), 1 + <M«2 - 1))(U2 - If 

Hence, 

wi d l(x, Ul ,u 2 ) = Ci d \x)( Ul l)y(x) k + d + C {d l h {x){u 2 - l)y(x) k + h + d 
+ 0((k + h)y(x) k + d (\ Ul - 1| 2 + \u 2 - 1| 2 )), 

where we can neglect the mixed derivatives as either (ui — l) 2 or (u 2 — l) 2 will determine the 
dominant part. For the second part we again use a Taylor polynomial, usin the fact that \x l \ < 
p < 1 and \u l r — 1| < i|u r — 1| for i > 2, r — 1, 2, hence the result follows immediately. □ 

Note that the terms ui — 1 and w 2 — 1 are asymptotically proportional: ^2. = e ^'~ 1 <~ and 
that y( x 2 ) k + h+d is exponentially smaller than y(x 2 ) k+d as h — £\/n- 
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Lemma 5.4. There exist constants such that Wjfl — w)lI{x,Ui,U2) is given by 
„„(<*) 



S d ) , 



w { d) h y(x) k 



k,h ,,. (d) 

W^-^^gf + o(K-i| + |n 2 -i| 

for Hi £ Ht,«2 £ Z'k+h a>nd x G A c , where fjf is given by: 

fc-i 



/fc (a;,Mi,M 2 ) = u^(x,mi,u 2 ) E 



S ; , h (x,wi,w 2 )y(a;) 



(37) 



Proof. We can argue similarly as in the proof of Lemma 14.51 and derive the recursive description 



(d) (d) , . n k,h 

W k+l,h = V W kM 1 + 



(d) 



and equivalently 



(d) 
w fc+i,h 



.(d) ~ 
J k,h 

~W 

k.h , 



w 



(d) 

k.h 



W 



Id) > 
J k,h 

Id) 



k,h , 



Further we get 

yk+l 



(d) 

u k+l,h 



(d) 
U k,h 



yW . 1 ,(T~\ k+1 1 



, w ki , 



Solving the recurrence leads to 



y 



i 



w 



(d) 

k.h 



,(d) 



^ -.,(<*)„„(<*) 2 1-y 



fc-i 



<?(E*w 



' =0(^0,^) 



u y 



J l,h 



Id) 

l-« ^l,h 
=Q(£<) 



(d) 



fc-l 



1 — w, 



(d) \ - 
Q.h 



V+i 



id) 



,(d)(d) 



w o,fe i-r , nl 2(d) , /0/ (d) i - L fc 



V 



:fjf' (x,«i,u 2 ) 



/ 



Observe that 



= «) + («" l)^(5d_i, (x, «),..., yJV. »*)) - 1,(3!) 

= w {d) + (u - l)xZ{S d . x , y { h d) (x,v),..., y {d \x\ «<)) 

= Cl d \x)(u 2 - l)y(x) h+d + (ui - l)y{x) d = 0{\u x - 1| + \u 2 - 1|) 



□ 



In the following, we denote by U := (ui — l)y{x) d and W := w^L (x, u^)- Note that wq^ 
U + W. By LemmaPwe obtain for wfl (note that Cf^ h {x) = C id \x)(l + L 1 )) 



l i(x, Ul ,u 2 ) = Cf\x){ ui - l)y(xY +d + C^ h (x)(u 2 - l)y(x) 



t+d 



\t+h+d 



+ 0((£ + h)y(xY(\ Ul -l\ 2 + \u 2 -l\ 2 )) 
= y{xf (Cfu + C [d \x){u 2 - l)y(x) h+d 

+0((£ + h)(\u 1 -l\ 2 + \u 2 -l\ 2 ))) 
= y{xf{C {d \x)U + W)(l + 0(h(\ Ul - 1| + K - 1|)) 
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We use the representation (jzTT)) for Cg(x), which we already used in the proof of Lemma [4.51 
and omit all error terms, to obtain by telescoping 



fc-i c, 



(*>(■£?) 



Z+d 



^(C ; (rf) ( a; )t/ + ^)(C^ 1 (x)C/ + ^) 
-1 /^(rf) 

(d) / \ V i+1 



fe_1 (CffifaOtf + WO - (C7 ; (rf) (a;)C/ + W) 
^ (C^(z)[/ + W0(O^ + W0 



As we know from (IT51) , Cg = ^^r-^ = 0(1) near u = 1 (analytic), hence 



Using 



Mr, ^ f 1 (U + W) 

ci d \x)U + W / 



(ui - 1) ~ 4=' 



1 - - b 
and w^^x, Ui, u^j ~ W / + 17, we can derive 



w^ d) h y(x) k 



Cp d U+W 2 1-V(x) 



" v^-((cp^+4 d) (^«))^|a-e- Kb v^i) 

~~ ^Se^^ _ ((Cpt$+wj?>(x,v))^(tBi&(%by/=p8)) 

and with the expansion (HO of toj^a.Ua) with «2 = and ft, = S,y/n, given by Theorem [SJ we 
can derive the expansion given in Theorem |6] □ 

Proof of Theorem [H The characteristic function of the two dimensional distribution is given by 

yn 

1 [(d), 4 4, dx 



1+ 2^/X^ e ~' e ~W (38) 
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We use the same contour as in the one dimensional case. With the same arguments, only 
integration over 7 contributes to the result, hence the representation (|32|) of leads to: 



»$,„(*!>*») = 1 + 



1 + 2 log 71 



Cp d 
b^Tp 



i[h + 



_ e< i 5i)v ^ ) _ i£2 cpl sinh ( i ?bv ^ )) 



,(-jKbV-P s ) 



l-ilog 2 n -y/— S( 



1 ^Tp l tl + 



„(- Aei>^=p^) 



sinh ps)) 



(sinh (f 6 1 /=ps)) 



where Vv? (/i; ^2) is the characteristic function of the random variable ^ P - b u (^/§ K ' i^f^) ) ■ ^ 

6. Tightness 

We must show the estimate ([6]) in Theorem |3] The fourth moment in (|6]) can be obtained by 
by applying the operator \u-§^) and setting u = 1 afterwards. Hence, using the transfer lemma 
of Flajolet and Odlyzko [3U] it turns out that it suffices to show that 



d_ 

On 



d 3 d 4 \ 

7 — + 6 ^3+^l) yr,H(x,U,U-) 



du 2 



O 



(39) 



uniformly for x £ A and h > 1 (see jTTl pp.2046] for the detailed argument) 
Set 



{dm 1 \ 

T k (.x) 



d 3 y k {x,u) 



and 



<9%^ (a;,M, ^) 



The left-hand side of (|3"9")l is a linear combination of Jjfj^\x) for j = 1,2,3,4. Therefore we 
need bound for those quantities. We will derive upper bounds for all j since this more general 
result is easier to achieve. We start with an auxiliary result. 

Lemma 6.1. Let j be a positive integer. Under the assumption that for all i < j the bound 

r yjf^^ l \x) — O (\x/p\ k ) holds uniformly for \x\ < p, we have (-§-) J Sjj. = O {L k ) for some 

L J u— 1 

positive constant L < 1. 

Proof. By Faa di Bruno's formula we have 

d ^ 



du k 



T.7 



u=l i>2 



=1 

(d)[ vi - 



A=l 



A! 



By our assumption we have 7^, rf '" lyi+ +I/m '(a; l ,l) = 0(\x l /p\ k ). The product is essentially a 
derivative of order j = J2\v\ of u l and can therefore be estimate by 0(i^). So the whole 
expression is bounded by a constant times ~^2 i>2 V~ 1 x lk jp k = O ((\x 2 \/p) k i) = 0((p + e) k ). 
Hence we can choose L — p + e to get the desired bound. □ 

Exactly the same line of arguments yield the analogous result for two levels: 
Lemma 6.2. Let j be a positive integer and set 



(40) 



i>2 
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Under the assumption that for all i < j the bound — O (|x/ j o| fc ) holds uniformly for 



\x\ <pwe have {-§^) J 



u=l 



O (L fe ) for some positive constant L < 1. 



With the auxiliary lemmas we can easily get bounds for ~fu®^(x) and l^^Kx). 
Lemma 6.3. We have 



MM, 



Ik ( X ) 



O (1) uniformly for x € A, 

(|x/p| fc ) uniformly for \x\ < p 



(41) 



and /or £ > 1 



Ik \ x ) 



Ofminffc^ 1 ,^^ 
0(|x/p| fe ) 



(42) 



uniformly for x € A, 
uniformly for \x\ < p 

Proof. The estimate (j4Tj) essentially follows from Lemma FOI We know 7^ ^\x) = C^ d \x)y(x) k+d (l- 
O (L k )) = C (1) with some < L < 1 and |y(a;)| < 1 and this is sufficient to show the first part 
of CH]). 

If \x\ < p we can exploit the convexity of y(x) on the positive real line to get |y(ic)| < \x/p\. 
This implies l^ix) = O (\x/p\ k+d ), an even better bound than stated in the assertion. 

Now we are left with the induction step. Again we use Faa di Bruno's formula and the fact 
that u>j^(x, 1) = 1) = and obtain 



lk ( x ) 



9 {d)< \ 
— W K k '(X,U) 

ou 



n 



E A,!---.\,! 



= v{x) 



^ exp (x, u) + Ej^ (x, u) 



u=l 



d_ 

du 



y( x ) 



E 



n 



r V7»*) + r&"(*)r 



+ v(*)(7jS?i M (*) + rn* J ) ( 4 3) 



where T 



fc-i 



A,!---A,! . . 



Consider the case |x| < p. The product comprises only terms which essentially have the form 

l ( k-i ] ( x ) + T k d -i ] ( x ) with 3 < L Thus b y the induction hypothesis, iflf'ix) = 0(\x/p\ j ). 
Therefore the assumption of Lemma 16.11 is satisfied and the terms as a whole are bounded by 
C ■ \x/p\i. Since Ylj=\j^j — & we get 

7? M (*) = *(*)(7ffi M (*) + rSS q + o(l*/p| < )- 

So we finally get the desired estimate by induction on k and Lemma 16.11 starting with 

{dm = ixZ d ^ 1 (y(x),y(x 2 ), . . . ,y( x<i ~ 1 )) if ^ = 1, 
I else. 



7o 



(44) 



Now let us turn to general x G A. Like before we focus first on the terms of the product of 
43]) . Again the induction hypothesis guarantees that the assumption of Lemma 16.11 is satisfied 



and so rl_j (x) is exponentially small. Furthermore, the induction hypothesis implies jjf_i(x) = 
O (min Since j ( ^ ] (x) = O (1) this implies 
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Set 



(l-|y(z)|)£A^ 

Note that X^=iC/ — = ^ — 2j=i ^j- Since the term correspondung to \g — in Faa di Bruno 
formula is the very last term in (|4"5)l . we must have Jj^Zj — ^ and thus A < k l ~ 2 . Moreover, 



we have 



-l 



- 2 )Aj = £- h - k o < £- 3 
since J^i=2 < ^ implies fci > and, in particular, Y^j=2 -\j = 2 implies fci = I. Therefore 



B < 



■ij- 

k i-3 



We want to show that 

5 < - (46) 

Set Aj = and = fc^ 2 /(l - |y(a;)|). Note that Bj < Aj is equivalent to 1/(1 - \y(x)\) < k. 
Therefore the term B appears in our upper bound if and only if x is such that 1/(1— \y(x) |) < k. 
But this implies that B < prn^ji as desired, because the desired bound is equivalent to 

1 < fc(-2+2£'=i\,->i)/(-l+E'=2\f) = fc i+« 

1 - 



where a = (X) J= i — 1J / 1 + Sj=2 J > an< ^ ncncc the desired bound (|4^|) is weaker 
than 1/(1 - |j/(ar)|) < k. 

Now let aj; := Ji? (x). We have shown so far that 

a k = y(x)a k -i + y{x)A k with A k = O [ min k l ~ 2 , 

V V i-ls/wl 

and we know that ao is given by (j44[) . Solving this recurrence relation gives 



a/c = y(a;) fc a + C 



2/0) 



1 - y(x) 



1 - y(x) 



min ( fc £ 2 



fc 



£-3 



Since 



y(x) \ < and aoy(x) k — O (y(x) k+d ^ = O (1) we get the desired bound for a k and 



, i-p(gr 

l-y(x) i 

the proof is complete. □ 
Lemma 6.4. We /iaue 



and /or t > 1 



^WMm _ / W uniformly for x E A, 

7fc,fc W | C(|a;/p| fc ) uniformly for \x\ < p, [ ' 

7 ( d )M( x ) = / ( min i-|t/(1)| )) uniformly for x E A, 

[ O (|x/p| fc ) uniformly for \x\ < p 

Proof. As the bounds are precisely the same as in the previous lemma, the induction step works 
in an analogous way, using Lemma 16.21 instead of Lemma 16.11 Thus we only have to show the 
initial step of the induction, Eq. (|47|) . 

We can use a similar reasoning as in the proof of [171 Lemma 7]. Indeed, by applying the 
operator [ Jjj ■ ] _ x to © we obtain the recurrence relation 
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with initial value 7q < ?^ (se) = xZd-i(y(x), y(x 2 ), . . . , y{x d ~ 1 )) — 7/1(2:). Induction on k gives the 

representation 7^ [1] (x) = 7 £ d)[1] (x) - 7^ (x) and using 7 ^ )[1] (x) = + (L k )) 

from Lemma l4.1l we obtain 



7i d)[11 (*) = O (sup |y(x) fe+d (l - y(x) h )\ + L k ^j = O (- 



Since the last term is bounded, the proof is complete. □ 
Now, applying Lemma \6. 41 to (|39|) proves tightness and Theorem [3] after all. 

7. The Joint Distribution of Two Degrees 

We want to gain knowledge on the correlation between two different degrees d\ ,6,% in a certain 
level k — Ky/n. 

7.1. The covariance Cov(Xri il \k), xi d2 \k)). The covariance of two random variables X and 
Y is given by 

Cov(X,Y) = E{XY) - E(X)E(Y). 

In this section, we will prove the result on the covariance function of the two random variables 
X^ 1 \k) and X^f 2 \k), counting the vertices of degree d\ and e?2, respectively, on level k, given in 
Proposition 12.11 

To compute E (xi dl \k) ■ X^ 2 \k))\ we need to determine 



— [x r - 

Vn 



dudv 

while E(xi dl \k)) and E(X^ 2) (fc)) are given by 



d 2 

y k (x,u,v) 



We use the notations 



-—\ xU ]l^k 1 \ x ) an d — \ xn Wk 2 \ x ) i respectively. 



(di)/ \ 9 , u 2 \ d 

Tk {x,u,v) = —y k (x,u,v), 7^ '(x,u,v) = — y k {x,u,v), 

as well as (recall [TT] 

f fc <*'*5(z,u,v) = Jf-y k (x,u,v) andffc (dl ' d2) W =7fc (dl ' d2) (^l,l)- 
out; 

Lemma 7.1. There exist constants e and 9 such that for z 6 A(rj,6) 

~^\x) = C<*> (a: ) • C(*)(x)tf(a:)^+* + O(L')), 

where C^ dl \x) and C^ d2 \x) are given in Lemma \4-1\ 

Proof. We use the recursive representation (|20[) for -y( dl ) (a;, u, v) with the additional variable v. 
This gives 

(rfl)/ \ (d) 1 \V (di)/ i i i\ i—1 

% + l{x,u,v) = 2/fe+i(a;,M,u)2^7fc ; (a;,u,u)ii . 
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Derivating with respect to v gives 

^i d2 \x,u,v) = ^{x^v^l^^yyy- 1 

»>1 

~ (di,d2) / i i „.i\-A— 1 
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y fc+ i(a;,w,u)^i7^ 1 ' <i2) ( 



x ,u ,v )u v 



= y k+1 (x,u,v) ( ^^V^XK" 1 ^'^'"V 

,i>l / \i>l 

£>1 



+ 2/fc+i(z,u, v)2_^i% J (x,u,v)u v 



with 7Q dl ' d2 \x) = 0. Setting u = v = 1 we obtain 



lk : +l ^-\x)=y{x) 



K i>l J \i>l J i>l 

= Vix) {{li dl \x)+T^\x)){ 1 ^\x)+T^\x))+t^ , 
where we use the notations rl rfl ^(x) = 7fe dl ' > ( x ') an d r^ 2 ^(x) = lk' 2 \ xl ) as i* 1 the proof of 

i>2 i>2 

Lemma FTT1 and T^ ud *\x) = £7^ (a:*). Solving the recurrence, we get 



i>2 



il dud2) (*) = X>(^ ((7i dl) (») + r^ l) (x))( 7 ^ ) (x) + ri^(x)) + ff (48) 

From Corollary E3] we know that Tf l] {x) = 0(L e ) and r^ 2) (a;) = 0(L e ) in 9 for u = v = 1. 
Together with Equation (fT6| we have 



fc-i 



l+d 2 



0(L e ))+0(L l 



and the result follows. 



□ 



To extract coefficients we will use Cauchy's formula, 
where 8 is the truncated contour 8 = 8i U 82 U £3 U £4 given by 



x ) > P' 1 

n 



6 2 = < x = p ( 1 



p < a < p + rj log 2 n/n 





7T 












~ 2 



(49) 



$3 = < x = a 



p < a < p + rj log 2 n/n 



and 84 being a circular arc closing the contour, cf. Figure [2j 

It can be shown that the contribution of the circular arc £4 is exponentially small and thus 
asymptotically negligible. Near p, more precicely for z — p(l + — ) G 81 U 82 U #3 and = K-y/n, we 
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(5-contour 



FIGURE 2. The integration contour S 



have 



y( X fl+d2 =l + Q 

1 - y(x) ~ b^fp 



1 + 



y{xf ~ exp(-K&V^ps) ( 1 + O 
C^\x) ~ C^p* + , C^\x) ~ C d2 / 2 + O (|^|) . 

Hence, the expected value [x n ]^ dl ' d2 \x) is given by 



1 



2iri Jg by/p\/—s n 



as X^=o = 1-^(2) • Integrals of the shape f s (— s)' i e s ds can be easily transformed 

into Hankel's representation of the Gamma- function and together with y n ~ ^£^-"77,— V 2 ( c f 
Equation ([3])) and a = nb^fp we obtain 

E (X^{k) ■ xM(k)j) = C dl C d2 p d ^^-i 



2 i2 

-K b p 



From this the representation of the covariance given in Proposition 12.11 easily follows. 

7.2. The correlation coefficient. To obtain more information on the correlation of two degrees 
d\ and d 2 on the same level fc = n^/n, we compute the correlation coefficient, given by 

o.(**) W ,*s« W )) . ^m^m . 

^Var(x( dl) (fc)) V /Var(^ 2) (fc)) 
For the computation, it remains to compute the variance given by 

Var(X^)(fc)) =E((X^)(fc)) 2 ) - (E(X^(k))) 2 . 
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Figure 3. The covariance ar(iJ„) 



We need to determine E ((X^ 1 \k)) 2 ^ , which can be done very similarly to the previous part. 



9 ( 9 ( U 



u=l 



Proposition 7.2. The Variance Yax(X^ (k)) of the random variable Xn dl \k) counting vertices 
of degree d\, with d\ fixed, at level k = n^/n in a random Poly a tree of size n is asymptotically 
given by 



Yav(X^>(k)) = C dl C d2 p^n 
as n tends to infinity. 



b 2 p 



-K 2 b 2 p 



K e 



0(Vn), 



(50) 



We proceed analogously to the computation of the variance, and obtain the following auxiliary 
result. 

Lemma 7.3. There exist constants e and 9 such that for z £ A(n,0) 
^P]>(s) = (C<«« (:b) ^ 

1- 2/0*0 

where C {dl \x) is given in Lemma \4-1\ 

Proof. The proof of this lemma is analogous to the proof of Lemma EU derivating recurrence (|20l) 
a second time. The additional summand C^ dl \x)y(x) k+dl origins in derivating twice with respect 
to the same variable u. □ 

Note that the additional summand C^ dl \x)y{x) k+dl in Lemma 17.31 where ji dl ^\x) and 
li dl ' d2 \x) differ from each other, is equal to the expexted value ¥,(x^ 1 \k)\ when extract- 
ing coefficients -^-[x n ]C^ dl \x)y{x) k+dl . As this is of order \fn, while the coefficient of the other 
terms will be of order n, this term is negligible, and we obtain 

2 / _«V£ 



E 



({X^\k)f)=C dl C d2 p 2d ^ 



b 2 p 



2 i2 

-K b p 



(51) 



by using Cauchy's formula and the integration contour S given in (|49[) . Applying the known 
estimate for Epff, (k)) we obtain the representation given in Proposition 17. 2\ and with Propo- 
sition 12.11 the result given in Theorem 0] follows immediately. 
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